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ABSTRACT 

Due to stellar rotation, the observed radial velocity of a star varies during the transit of a planet 
across its surface, a phenomenon known as the Rossiter-McLaughlin (RM) effect. The amplitude of the 
RM effect is related to the radius of the planet which, because of differential absorption in the planetary 
atmosphere, depends on wavelength. Therefore, the wavelength-dependent RM effect can be used to 
probe the planetary atmosphere. We measure for the first time the RM effect of the Earth transiting the Sun 
using a lunar eclipse observed with the ESO HARPS spectrograph. We analyze the observed RM effect at 
different wavelengths to obtain the transmission spectrum of the Earth’s atmosphere after the correction of 
the solar limb-darkening and the convective blueshift. The ozone Chappuis band absorption as well as the 
Rayleigh scattering features are clearly detectable with this technique. Our observation demonstrates that 
the RM effect can be an effective technique for exoplanet atmosphere characterization. Its particular asset 
is that photometric reference stars are not required, circumventing the principal challenge for transmission 
spectroscopy studies of exoplanet atmospheres using large ground-based telescopes. 

Subject headings: stars: rotation — eclipses — Earth — planets and satellites: atmospheres — techniques: radial 
velocities 


1. Introduction 

Among the nearly 2000 exoplanets discovered so 
far, more than half are transiting systems. Due to the 
rotation of the host star, the observed stellar radial ve¬ 
locity (RV) is expected to change as the planet transits 
different parts of the rotating stellar surface. This is 
called the Rossiter-McLaughlin (RM) effect and was 
initially measured for eclipsing binary stars dRossitei 
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1924 : McLaughlinl 1924 ). After the first observatio n 


of the exoplanetary RM effect by Oueloz et al. (2000), 
there are now more than 80 exoplanet systems with 
this effect observed, e. g. Albrecht et al. ( 2012 ). This 
effect is normally used to measure the projected angle 
between the stellar spin and the planetary orbit. 

Snellen ( 20041) was the first to propose using the 


wavelength-dependent RM effect to probe the atmo¬ 
spheres of transiting exoplanets. This method ex¬ 
ploits the fact that the planetary effective radius is 
wavelength-dependent due to differential atmospheric 
absorptions. Dreizler et alj ( 2009!) theoretically mod¬ 
eled the RM effect for the Na absorption of giant ex¬ 
oplanets. The advantage of the RM method compared 
to the traditional spectrophotometry is that it does not 
demand the use of photometric reference stars. The 
RM method promises to become a powerful technique 
for future transmission spectrum measurements, es¬ 
pecially as the next generation of very large ground- 
based telescopes are likely to have a small field-of- 
view, making it difficult to employ suitable reference 
stars. 
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With more and more terrestrial exoplanets being 
discovered, the characterization of their atmospheres 
will become a major goal for exoplanet research in 
which the RM method can play an important role. Al¬ 
though it is currently difficult to observe the RM effect 
of these terrestrial exoplanets due to instrumental limi¬ 
tations, a lunar eclipse provides us with an opportunity 
to observe the Earth transiting the Sun and so explore 
the effectiveness of the RM method. The concept of re¬ 
garding the Earth as an exo-Earth using lunar eclipses 


has_ been app lied in the past by [Palle et al 
VrdaEMadjtffetal. ( 201(1 ) 


( 2009 ): 
U golnikov et all ( 2013 ): 


Arnold et al.1 ( 2014 ): Yan et al. d2015b ). These studies 


obtained the transmission spectrum directly from the 
ratio spectrum before and during a lunar eclipse. Here 
we present for the first time the use of the RM effect 
of an Earth transit during a lunar eclipse to retrieve the 
transmission spectrum. 


2. Observations and RM effect measurement 

We observed the 15-April-2014 lunar eclipse with 
the High Accuracy Radial velocity Planet Searcher 
(HARPS) mounted on the ESOLa Silla 3.6m telescope 
( Mayor et al.ll2003) . A consecutive sequence of obser¬ 
vations over one entire night covered all of the eclipse 
stages, i.e. the penumbral and umbral eclipse and out 
of eclipse (called hereafter the bright Moon). The clas¬ 
sical fiber spectroscopy mode was used with the fiber 
located at the center of the Tycho crater (see Fig.[T|for 
the trajectory). The telescope used non-sidereal track¬ 
ing and the tracking velocities were updated every few 
minutes. In addition, manual guiding on crater struc¬ 
tures was performed. This worked reasonably well and 
limited the pointing coordinate drift to < 1.5 arcmin. 
We used varying exposure times during the eclipse be¬ 
cause the lunar surface brightness changed dramati¬ 
cally. In total, 382 lunar spectra were obtained (Table 

m. 

Data reduction was performed using the HARPS 
pipeline. Each exposure frame comprises 72 spectral 
orders that cover the wavelength range from 378 nm to 
691 nm with a spectral resolution of A/AA ~ 115,000. 
Radial velocities are measured, for each spectral order 
as well as for the overall spectral range, using a cross¬ 
correlation with a G2 stellar spectral template. The 
overall measured and theoretical RV curves are shown 
in Fig. [2] The theoretical RV is the combination of 
two components: the motion of the Sun with respect 
to Tycho and the motion of Tycho with respect to the 


observer. These velocities are calculated with the JPL 
Horizon Ephemeri'H, which considers the orbital mo¬ 
tions, rotations of the Earth and the Moon, and the light 
travel time corrections. 


There is an offset between the theoretical RV and 
the measured RV, which probably originates from the 
RV zero-point of the spe ctral template used by the 
HARPS pipeline ( Molaro et al.ll2013 ). We correct this 
offset by using the data points taken during the bright 
Moon as a baseline. Fig. [2] shows the measured RV 
with the motion RV corrected and the baseline offset 
subtracted. This is the final RV curve of the RM ef¬ 
fect for the Earth transiting the Sun. At the bright 
Moon stage, the RV is essentially corrected to zero, but 
there is still a small slope with an amplitude of about 
4 m/s which is probably due to the instrumental drift 
and the non-perfect telescope guiding when observing 
the Moon. Since the RM amplitude is on the order of 
2 km/s, this small residual slope does not significantly 
affect our analysis. 


When the Moon enters the 1st penumbra, the ob¬ 
served RV becomes negative (blueshift) as the Earth 
begins to transit the redshifted rotating part of the solar 
disk. At the umbral stage, the RV gradually changes 
from a blueshift to a redshift. When the Moon enters 
the 2nd penumbral stage, the RV is redshifted since 
the Earth obscures mainly the blueshifted solar region. 
The RV gradually decreases during this stage as the 
Moon moves out of the penumbral shadow. 

The umbral RV has a relatively complex structure 
and its details are determined by a combination of both 
the refracted part of the solar disk and the properties of 
the Earth’s atmosphere refracting the sunlight. In this 
work, the umbral part is not studied since we focus 
only on the penumbral parts for the RM effect. 


3. Retrieving the transmission spectrum 
3.1. Modelling the RM effect 


To retrieve the wavelength-dependent RM effect 
due to the Earth’s atmospheric absorption, we firstly 
establish an RM model to take into account the 
wavelength-dependent solar parameters. 


Ou r mo del follows the m ethod described in 'Gaudi & Winnl 
(2007) and Haswell ( 2010h . We divide the solar disk 
into elements with a size of 0.01R 0 x 0.01R 0 and the 
radial velocity due to solar rotation for each element 


-Http://ssd.jpl.nasa.gov/?ephemerides 
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Fig. 1.— Trajectory of the observed position (Tycho crater). Our observation began at the bright Moon stage and 
lasted until the 2nd penumbra. The figure is reproduced from the NASA lunar eclipse page. 


Table 1: Details of the observed spectra. 


Eclipse stage 

UT time span 

Number of observed spectra 

Exposure time (second) 

Airmass range 

Bright Moon 

00:30 - 05:32 

294 

30 

1.98 ~ 1.09 

1 st Penumbra 

05:32 - 06:45 

55 

30,100,500 

1.09 ~ 1.23 

Umbra 

06:45 - 08:47 

13 

500 

1.23 ~ 2.01 

2nd Penumbra 

08:47 - 09:35 

20 

500 ~ 60 

2.01 ~ 2.90 



Fig. 2.— The observed RV (red line) and the theoretical RV of the motions between the Earth, the Moon and the Sun 
(black line). The observed RV after the corrections of the motions and the baseline offset (green line), which is the 
RM curve of the Earth transiting the Sun. Positive values represent redshift. 
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Fig. 3.— The observed RM effect curve (green points) together with the modelled RV curve (red line), which is the 
combination of the white-light RM model (black line) and the convective blueshift model (blue line). See the text for 
the details of the parameters used in the model. Two regions on the RV curve which are symmetric with respect to the 
mid-eclipse (labelled as PM1 and PM2) are chosen to calculate the RM amplitude. PM1 contains 4 observed points 
with a time span of 496 s and PM2 contains one observed point with a time span of 500 s. The lower panel shows the 
residual between the observed and the modelled RV. 
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is: 


v(x, y) = oj x R 0 sin I s (1) 


where u> is the rotation angular velocity, I s is the incli¬ 
nation of the solar spin axis towards Tycho, x and y are 
in the sky-plane coordinate system of which the origin 
is at the projected solar center and the y-axis is along 
the projected solar rotation axis. The coordinate values 
are in units of the solar radius R G . The RV anomaly 
A V(t) due to the RM effect is then calculated by inte¬ 
grating the intensity-weighted RVs of the visible solar 
disk: 


A V(t) = 


f f v(x, y ) I(x, y) d.r dy 
f f I(x, y) dv d y 


( 2 ) 


The RM curve for white-light (i.e. the entire wave¬ 
length range of the HARPS spectrum) is generated 
from this model (the black line in Fig. [3]). The pa¬ 
rameters are further described as follows: 


1. the coordinates of the Earth’s center and the so¬ 
lar center as seen from Tycho are generated us¬ 
ing the Horizon Ephemeris. We use the Earth 
radius of R® = 6378 km and the solar radius of 
R 0 = 6.955 x 10 5 km. These data are used to 
determine the visible part of the solar disk at a 
given time. 


2. the inclination I s is calculated using the solar 
north pole position angle given by Horizon. The 
mean value of I s is 95°.61 and changes by about 
0°.02 during the eclipse. 


3. the solar differential rotation used is (iCoxl2000l) : 


a> - A + Bsirrcp (deg/day), (3) 


where </> is the heliographic latitude and is calcu¬ 
lated with: 


sin <p — (1 - x 2 - y 2 ) 1 ^ 2 cos I s + y sin I s . (4) 


The coefficients adopted are A = 13.46, B = 
2.99 JPupree & Henzell 19721) . 


4. the quadratic limb-darkening coefficients (Claret 


2004) are used for the white-light model to 


calculate the intensit y I(x, y). W e adopt the 
coefficients used by Moiar o et at. (120131) . i.e. 
u a = 0.5524 and u * = 0.3637. 


However, we adopt the actual data from the Horizon 
Ephemeris and so the model curve in Fig. [3] is not the 
result of such a fit. The “projected spin-orbit” angle for 
this observation is -169°.7 calculated using the trajec¬ 
tory of the Earth as seen from Tycho. Because this 
trajectory is the combination of the Earth and lunar or¬ 
bits, the angle is not the real spin-orbit angle of the 
Sun-Earth system of 7°.25. 


3.2. Other effects 

3.2.1. Convective blueshift 


The stellar convective motion results in a blueshift 
of the lines which is called the convective blueshift 
(CB) ( Dravins 1982 ). As the CB radial component 
varies across the stellar disk, the observed stellar RV 
changes when a planet transits the stellar surface. 


Shpo rer & Brownl ([2011) discusses the CB during an 


exoplanet transit and its effect on the RM effect. We 
follow the CB model in Shporer & Brown ( 201 ll) and 
adopt a typical CB value of the Sun Vch ~ -300 m/s 
for the white-light RV calculation. The blueshift value 
for each element on the solar disk is pVcs, where 
p = (1 - x 2 - y 2 ) 1/2 . The RV change caused by CB 
(blue line in Fig. [3} is then calculated in a similar way 
to the RM model and added to the RV value caused by 
the RM effect to give the red line in Fig. [3] 


3.2.2. Atmospheric refraction 


The effects of atmospheric refraction are impor¬ 
tant for interpreting the transits of terrestrial planets 


(Garci a Mufioz et al ■ 2012: iBetremieux & Kalteneggei ' 


2014) and affect the RM effect. As the atmosphere re¬ 
fracts the light from the stellar disk region that should 
be obscured by the planet, the apparent effective ra¬ 
dius of the planet will appear smaller than without 
refraction. Thus in general, refraction would result 
in a smaller RM amplitude than without refraction. 
In our RM model, we do not include atmospheric re¬ 
fraction, which means the modelled RM amplitude 
should be larger than the actual observed RM ampli¬ 
tude. This is consistent with the result in Fig. [3] in 
which the modelled RV value without refraction (red 
line) is generally larger than the observed RV (green 
points). 


3.2.3. Line broadening mechanisms 


For an exoplanet transit, the projected angle be¬ 
tween the stellar spin and the planetary orbit is calcu¬ 
lated by fitting the observed RV with the RM model. 
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In the RM model described above, we assume the 
line profile is dominated by the stellar rotation and 






































use t h e first-moment a pproximation dGaudi & W inn 
20071: lOhta et al.ll BoQ5). This assumption gives us 
the RV anomaly expressed by Equation [2] However, 
other mechanisms like the microturbulence, macrotur¬ 
bulence and instrumental broadening also affect the 
line profile and can change the shape of the RV curve 
(Hirano et al. 2011). This may contribute to the dis¬ 
crepancy between the modelled RV and the observed 
RV in Fig. [3] 

3.3. The RM effect in different wavelength ranges 

The HARPS pipeline produces directly a RV value 
for each of the 72 spectral orders. Three orders do 
not have useful RV data due to the lack of stellar lines 
in the corresponding spectral mask. After correcting 
for the Earth-Moon-Sun motion and the baseline offset 
described in Section 2, we obtain 69 RM curves each 
representing a different wavelength range. 

To obtain the RM effect differences caused by the 
Earth’s atmosphere, the wavelength-dependent solar 
parameters need to be considered. The first aspect 
concerns the stellar convection blueshift. The actual 
RV caused by CB depends on where the line forms, 
e.g. the strong lines or low-excitation lines may have 
a small CB because they usually form high in the stel- 
lar atmosphere w here the granulation is weak (IDravins 
1982 : Gravll2009l) . As each of the HARPS spectral or¬ 
ders comprises stellar lines with different CB values, 
the Vcb of each spectral order varies. Therefore, un¬ 
like the white-light CB model for which we can as¬ 
sume Vcb ~ -300 m/s, the actual CB effect for a given 
HARPS order is difficult to model. Instead, we use 
two symmetric parts on the RV curve to cancel this 
effect utilizing the fact that the convection induced RV 
anomaly during transit is symmetric with respect to the 
mid-transit as shown by the blue line in Fig. [3] Two 
regions of the RM curve (labelled as PM1 and PM2 in 
Fig. [3]) are then chosen. The PM1 and PM2 are sym¬ 
metric regions of the eclipse, i.e. they are at the same 
distance from mid-eclipse. Fig. |4jt shows the RV val¬ 
ues of PM1 and PM2 and the RV difference between 
them (i.e. PM2 - PM1). We use this PM2 - PM1 value 
to represent the RM effect amplitude for each spectral 
order. From Fig. [4jr, it can be seen that there is a cor¬ 
related systematic difference between PM1 and PM2. 
We believe that this correlated pattern results from the 
measured effect of the CB depending on the type and 
number of spectral lines present in each spectral order 
but is well-cancelled by using PM2 - PM1. 


The second wavelength-dependent solar param¬ 
eter is the limb-darkening. Here we use the em¬ 
pirical power-law limb-darkening coefficients from 
Hest roffer & Magnan d 19981) . For each of the HARPS 
orders, we interpolate a limb-darkening coefficient at 
the corresponding wavelength and calculate a value 
of PM2 - PM1 with a fixed Earth radius (6378 km) 
using our RM model. The model values of PM2 - 
PM1 are shown in Fig. St). The final RM effect at dif¬ 
ferent wavelength ranges, after the correction of the 
limb-darkening, is presented in Fig. St. The limb- 
darkening correction made here is mainly for the con¬ 
tinuum, however, since the solar lines have different 
limb-darkening compared to th e adjacent continuum 
( Snellen 2004 : Yan et al. 2015ah . this could introduce 
some extra noise in the final RM amplitudes. 

3.4. Results and discussion 

The RV curve of the wavelength-dependent RM 
amplitudes in FigSt results from the different effec¬ 
tive radius at different wavelengths, and it can be re¬ 
garded as a mapping of the low resolution transmis¬ 
sion spectrum of the Earth’s atmosphere into effec¬ 
tive radius and hence radial velocity space. To in¬ 
terpret the observed atmospheric features, we build 
a transmission spectral model foll owing the method s 
of Kaltenegger & Traubl ( 2009i) and Yan et al. ( 2015b ). 
We calculate the overall transmission spectrum of the 
Earth atmosphere from 0 to 80 km altitude considering 
the ozone absorption and Rayleigh scattering extinc¬ 
tion (blue line in Fig. [4}:). This transmission spectrum 
is overlaid with the RM amplitude curve to compare 
their shapes. We emphasize here that the transmission 
spectrum model is not used to fit the observed RM am¬ 
plitude curve but to demonstrate the presence of the at¬ 
mospheric features mapped into it. By comparing the 
spectral shapes in Fig. [4}c, we interpret the RM am¬ 
plitude curve as follows: the RM amplitude is larger 
towards the blue due to Rayleigh scattering extinction 
that makes the atmosphere more opaque and the at¬ 
mospheric effective thickness larger at shorter wave¬ 
lengths. The broad peak around 600 nm results from 
the ozone Chappuis band absorption while the RM am¬ 
plitude towards the red becomes smaller since both the 
ozone absorption and the Rayleigh scattering become 
weaker, rendering the atmosphere more transparent. 

The RM model we have used for this work suffers 
from several defects and incompleteness. The first is 
observational in that there is a slight, but barely signif¬ 
icant, drift in the telescope guiding on the lunar sur- 
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Fig. 4.— (a) Top panel: the RV values of PM1 (green line) and PM2 (red line) for the 69 spectral orders. The black 
line is the value of PM2 - PM1. All three lines here are shifted to zero at the first point (380 nm) in order to compare 
their structures. At the bottom, the RV curves of two bright Moon spectra are plotted (their values are shifted to - 60 
m/s), and they are relatively flat compared to PM1 and PM2. The order at 658 nm has a low RV accuracy due to a low 
number of available solar lines, (b) Middle panel: model values of PM2 - PM1 considering the wavelength-dependent 
limb-darkening coefficients with a fixed Earth radius (dashed line). The black line is the observed PM2 - PM1 as in 
the top panel, (c) Bottom panel: the difference between the observed and modelled PM2 - PM1 values shown in (b). 
This is the final RM amplitude in different wavelength ranges after the limb-darkening correction. The HARPS orders’ 
central wavelengths are indicated by asterisks. The transmission spectral model of the Earth’s atmosphere (blue line) 
is plotted for comparison. 
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face. The second is that other line broadening mech¬ 
anisms like the stellar micro/macro-turbulence are not 
included in the RM model (cf. subsection l3.2.3l) . The 
third is more fundamental and its full solution is be¬ 
yond the scope of this letter. This is the determination 
of the detailed mapping between the RM amplitude 
and the atmospheric exctinction as a function of wave¬ 
length. This depends on the effects of atmospheric re¬ 
fraction in the exoplanetary atmosphere which will be 
influenced not only by the absorption coefficients of its 
gaseous constituents but also by the effects of screen¬ 
ing due to clouds and aerosols. 

The combination of these effects can explain the 
slight deviation of the modelled RV from the observed 
RV values. It can also explain why the modelled 
RM amplitudes for the limb darkening correction are 
larger than the observed RM amplitudes (shown in 
Fig. BJ)). However, this has a very limited effect on 
our retrieval of the final RM amplitude curve since the 
limb-darkening correction produces only a small RV 
variation with wavelength so that the shape of the fi¬ 
nal RM amplitude curve is little affected by the RM 
model. Future modelling work containing a proper 
treatment of refraction and stellar turbulence will allow 
us to combine the transmission spectral model with 
the RM model, and to fit the observed RM amplitude 
curve directly instead of just comparing its shape with 
the transmission spectrum as presented in the current 
work. 

4. Conclusion 

We have observed for the first time the Rossiter- 
McLaughlin effect of the Earth transiting the Sun using 
a lunar eclipse. The RM effect curve has been obtained 
using high accuracy RV observations and an RM effect 
model is built to analyze the observed result. Separate 
RM curves at different wavelengths are obtained from 
69 HARPS spectral orders. After the correction of the 
wavelength-dependent limb-darkening of the Sun and 
the convective blue-shift of the solar lines, we retrieve 
the wavelength-dependent RM amplitudes due to the 
transmission of the Earth’s atmosphere. The ozone 
Chappuis band and the Rayleigh scattering signatures 
are clearly detectable. 

The RM method can be used to detect broad fea¬ 
tures, such as the Rayleigh scattering, in exoplanet at¬ 
mospheres. The advantage is that no reference stars 
are needed in contrast to the requirements for the tra¬ 
ditional spectrophotometric method. Since the next- 


generation ground-based telescopes such as E-ELT 
will have a relatively small field-of-view, limiting the 
access to nearby reference stars, the RM method will 
provide a promising technique for the characterization 
of planetary atmospheres. Thus, in the future, this 
method can be applied to detect the atmospheres of 
terrestrial planets and particularly to search for the bio¬ 
signature gas ozone. 

Extending the observation to near-infrared (NIR) 
wavelengths will yield more absorption features. 
However, the ability of the RM method is limited 
by the number of suitable spectral lines of the parent 
star. The Sun, for example, does not have sufficient 
lines in the NIR. However, for M-type stars which are 
expected to be promising targets for exoplanet atmo¬ 
sphere studies, there are many stellar molecular lines 
in the NIR which can be exploited. 
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